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MOMENT REPRESENTATIONS OF THE EXCEPTIONAL 
Xi-LAGUERRE ORTHOGONAL POLYNOMIALS 


CONSTANZE LIAW AND JOHN OSBORN 

Abstract. Exceptional orthogonal Laguerre polynomials can be viewed as an ex¬ 
tension of the classical Laguerre polynomials per excluding polynomials of certain 
order (s) from being eigenfunctions for the corresponding exceptional differential op¬ 
erator. We are interested in the (so-called) Type I Xi-Laguerre polynomial sequence 
degpn = n and a > 0, where the constant polynomial is omitted. 

We derive two representations for the polynomials in terms of moments by using 
determinants. The first representation in terms of the canonical moments is rather 
cumbersome. We introduce adjusted moments and find a second, more elegant for¬ 
mula. We deduce a recursion formula for the moments and the adjusted ones. The 
adjusted moments are also expressed via a generating function. We observe a certain 
detachedness of the first two moments from the others. 


1. Introduction 


Exceptional orthogonal polynomials were originally discovered as exact solutions 
to certain models in quantum mechanics. Those models include the supersymmetric 
setting, the Fokker-Planck, as well as the Dirac equations. See [5, 11, 12, 21, 24] for the 
connections to physics and mathematical physics. In recent years, the held has enjoyed 
much further attention from both the mathematics and the physics community; see 
e.g. [1, 4, 6, 7, 9, 10, 13, 17, 18, 22, 23] and the references therein. 

We focus on the exceptional Ai-Laguerre polynomials from the perspective of the 
seminal paper by Gomez-Ullate-Kamran-Milson [6]. This polynomial sequence, de¬ 
noted by a > 0, is orthogonal on [0, oo) with respect to the Ai-Laguerre 

weight 


lT“(a:) = 


{x + 

The polynomials are complete in L^([0, cxd); 1T“) even though there is no degree 0 poly¬ 
nomial. Further, they are the eigenfunctions of the exceptional Xi-Laguerre differential 
expression 


(l.l) t[y] = -xy” + [(x + a + l)y' - y]. 

\x + a J 

The spectral analysis of the polynomial system and a rigorous dehnition of a self- 
adjoint operator corresponding to was presented in Atia-Littlejohn-Stewart [1]. 
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The functions y{x) = T"(x) satisfy the eigenvalue equation 

f[y\ = (n — l)y (0 < a; < cxo). 


Probably the most comprehensive study of three types of exceptional Laguerre poly¬ 
nomials is contained in Liaw-Littlejohn-Milson-Stewart [17]. We refer the reader to 
Duran [4] for an interesting relation with so-called exceptional Meixner polynomials. 
Literature reveals several representations of the Xi-Laguerre polynomial sequence: 
For example, in [6] as well as Ho-Sasaki [13] the hrst couple of polynomials are listed 


( 1 . 2 ) 

(1.3) 

(1.4) 


L“(a:) = X + a + 1 , 

L^ix) = x"^ — a{a + 2), 

Lg (x) = - [—-f (a -|- 3)a:^ -1- a{a + 3)x — -\- Aa-\- 3)] , 


Further, these polynomials can be expressed in terms of the classical Laguerre poly¬ 
nomials {pnjnGNo (here p_i = 0), see e.g. [17, Equation (3.2)]: 

L‘^{x) = -{x + a + l)Kz!(x) -f- (x + a)pl_ 2 {x) (n e M). 

Duran’s work [4] contains representations of general exceptional orthogonal polynomials 
using determinants of classical polynomials. 

In [6] a three-term recurrence for the exceptional Wi-Laguerre polynomials was found 

0 = (n 1) [(x -1- af{n + a) - a]L"+ 2 (^) 

-I- (n -I- Q;)[n -I- a)‘^{x — 2n — a — l) + 2a]L^_^^{x) 

-1- (n -|- a — l)[(x -|- aYiji -|- a -|- 1) — q;]L“(x). 

We are most interested in the way the Xi-Laguerre polynomial sequence was in¬ 
troduced in the work of Gomez-Ullate-Kamran-Milson [8]. Namely, the sequence of 
polynomials {nj(x)}“g where 

(1.5) Vi{x) = X -|- a -|- 1 and Vi{x) = {x + a)* for i >2, 

span the exceptional Wi-Laguerre polynomial flag. Via the Gram-Schmidt process this 
sequence produces the sequence 

We note that the classical orthogonal polynomials are obtained from the sequence 
{l,x,x^,...} also by applying Gram-Schmidt. This perspective is used to find the 
classical moment representation, see e.g. [2, Ghapter 1.3]. An adaption of these ideas 
leads us to deduce our representations. 

The proof of our recursion formula for the moments relies on an application of so- 
called symmetry factors for differential equations (see e.g. Gole [3, p. 66]) which was 
further developed from second to higher order differential equations by Littlejohn, see 
e.g. [20, 19]. In essence, every second order differential equation is symmetrizable, and 
the corresponding symmetry equation is solved by the weight (of orthogonality). In 
Krall’s [15] well-known classihcation theorem (see e.g. Krall-Littlejohn [16, Theorem 
l(ii)]) this theory was used to derive a moment equation for the classical moments. 
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Here we carry out those techniques to obtain a recursive dehnition for the (adjusted) 
moments 



Jo 

A formula for the (exceptional) moments 



is obtained. 

Remark. We warn the reader that we often drop the subscript or superscript a to 
simplify notation. 

1.1. Plan of the paper. In Section 2 we present a preliminary result: We characterize 
the subspace spanned by the hrst m exceptional Xi-Laguerre polynomials (Lemma 
2.1), which yields what we call the exceptional condition satished by all exceptional 
ATi-Laguerre polynomials (Lemma 2.2). 

A hrst representation (Theorem 3.1) of the exceptional Wi-Laguerre polynomials 
in terms of the exceptional moments is found in Section 3. The expression is rather 
cumbersome. 

In Section 4 we focus on recursion formulae for the adjusted (Theorem 4.1) and 
the exceptional moments (Theorem 4.4). The adjusted moments are also expressed 
explicitly in two different ways (Corollary 4.2 and Theorem 4.3). One of the proofs 
involves a generating function. Throughout Subection 4.1 we notice that the hrst two 
moments are diherent in nature than the others. 

A more elegant representation (Theorem 5.1) of the exceptional Xi-Laguerre poly¬ 
nomials in terms of the adjusted moments is the topic of Section 5. In Remark 5.2 we 
determine the normalization constant so as to yield precisely those polynomials in the 
literature [6] and [17]. At the very end, we verify the representation formula for n = 1 
and n = 2. 

The author considers Theorems 4.3 and 5.1 to be the main results of this paper. 

2. The exceptional condition and the polynomials L" 

We characterize the span of the hrst m exceptional Xi-Laguerre polynomials as those 
polynomials of degree less than or equal to m for which the exceptional condition 


p\—a) — p{—a) = 0 


( 2 . 1 ) 


holds. This fact is well-known to specialists, but it is often written slightly diherent. 
We include the proof for the convenience of the reader. We mention aside that it is 
not hard to see that the polynomials given by (1.2) through (1.4), of course, satisfy 
the exceptional condition. 

In order to formulate precisely and prove the above characterization result, we let 
Vm denote the set of polynomials of degp < m and dehne the span of the hrst m 
exceptional Tfi-Laguerre polynomials 


Crn ■= span{L" : n = 1,..., m}. 
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and 


Mm ■= {p G Vm ■ P satisfies (2.1)}. 


Lemma 2.1. The sets Cm = Mm for all m G M. 


Proof. To show Cm ^ Aim take p G Cm- Then £"[p] G Cm, because Cm is a subspace 
of Lf‘{piP) which is invariant under In particular is a polynomial in x. In 

virtue of the dehnition of some cancellation has to take place in the exceptional 
term 

(^) [(^ + « + iy-p]. 

So the term in square brackets must equal zero when evaluated at x = —a. This yields 
precisely the exceptional condition (2.1). And so we have Cm ^ Mm- 

The opposite containment follows from a dimension argument. Namely, we have 
dim Pm = m+1. This implies dim Aim = rn, due to the imposition of the one restriction 
(2.1). We also have dim£m = rn, because the set is spanned by m polynomials with 
mutually different orders (implying their linear independence). Finally, we recall that 
Cm and Mm are both subspaces of Vm, and Cm ^ Mm- n 


The exceptional condition (2.1) yields a condition on the coefficients of the polyno¬ 
mials To see this, we write 

n 

Ln{x) = 

k=0 

and compute 

n 

{L^y (x) = 

k=l 

Since L" G Cn, condition (2.1) applies and we have 0 = (L^)' {—a) — L"(— a). We 
conclude: 


Lemma 2.2. The coefficients of the Xi-exceptional Laguerre polynomials L^{x) = 
YJk=o(^nkX^ obey 

n 

(2.2) -CnoA^ Cnk [k{-a)^~^ - (-a)^] = 0. 

k=l 


3. First Representation of L" 

Fix n E N. Recall that = X]fc=o ^nkX^- We determine Cnk for /c = 0,1, ..., n by 
means of a linear system of n -I- 1 equations Ac = b, where 


CnO 


■ 0 ■ 

Cnl 

_ ^nn . 

G and b : = 

0 

. . 


c : = 


G 
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with normalizing constant and A is given below. In this section, we do not £x 
Kn, but merely assume 7 ^ 0. In Remark 5.2 below we determine the normalization 
constant such that we obtain exactly the polynomials described in [ 6 ] and [17]. 


Theorem 3.1. The exceptional Xi-Laguerre polynomials admit the representation 

L^ix) = -7 > (det AAx^ = - -- -I 2 n (n e M), 

^ detA^^ ^ det A 1 x ... ^ 

k=0 

where 


A = 


Pi + (o + l)/io 


E 

m=0 


m 




° — f—a')^ 


l(-a) 
yU .2 + + 1 )/^! 


n —a 


in —1 


-a 


hn+1 + (o + l)Pr; 


2—m 


E 

771=0 


m 


/^m+l ^ 


2—m 


E 

m=0 


m 


/^m+n ^ 


2—m 



the exceptional moments are given by /i^ = /i^ = x^W°'(x)dx, and where the matrix 

Ak is obtained from A by replacing the (k + l)-st column with the vector b; as is done 
in Cramer’s rule. (In Section 4 below we find a recursion formula for the moments 
h'k-) 

Remark. There is, indeed, no polynomial of order zero. 

Remark. The matrix A is invertible, since the exceptional Xi-Laguerre polynomials are 
determined uniquely by exactly those conditions. Indeed, Lemma 2.1 ensures that the 
polynomial X]fc=o ^nkX^ belongs to the vector space and the other conditions given 
by the rows of the linear system simply require that it is orthogonal to the subspace 
Cn-i- Since dim \ Cn-i) = 1, the polynomial L" is uniquely (up to choosing the 
normalizing constant 0 ) dehned by these conditions. 

Proof. One condition on the coefficients Cnk was given by (2.2), which we rewrite in 
terms of the dot product 

[ —1 1(—a)° — (— 2(—aY — (—afi ... n(—a)"'~^ — (—a)” ] -0 = 0. 

We use the row vector in square brackets as the hrst row of the matrix A. 

The other n rows are obtained from orthogonality conditions: Recall that the poly¬ 
nomials are obtained from the sequence Vi given by (1.5) via the Gram-Schmidt 
process. In particular, we have 


(L“ , Vk)ivc = Kndnk for /c = 1 ,..., n. 
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where 6nk is the Kronecker delta. We now hnd explicit expressions for these conditions 
for the different values of A; = 1 ,..., n. 

For fc = 1, we have 

Kn^nl — {Ln i 

CnkX^ , X + a + 1 

\ k=Q I 

poo ! n n \ 

/ ^ + (a + 1) ^ CnkX^ ) hF“(a;) dx 

\k=0 k=0 J 

n 

Cnfc(/ifc+l + (a + l)/ifc). 




k=0 


The second row of the matrix A is determined by the first factor of the dot product 
[ /il + (a + l)po /i2 + (« + l)Ail • • • z^n+l + (« + l)/in ] ■ C = iFn^nl- 
When n > 2 we consider 2 < s < n and use the binomial theorem, 


Vg = {x + 


m=0 


m 


X™ a*-™. 


We compute 


KJr^s = {K , X,) 


W° 


-^0 \fc=0 / \m=0 k / / 

«C)0 ^ 

I E 

-^0 \m=0 




/c=0 


y^m=0 




k=0 


\m=0 


m 


m 


f^m+k 


Rewriting the summation as a dot product as before, we find that the (s + l)-st row 
of the matrix A is determined by 


y~! ( m ) ^ ( m ) ‘ ^ ^ ^ “ 

m=0 ^ ^ 7n=0 ^ ^ 771=0 

We solve the linear system Ac = b using Cramer’s rule 

Cnk = (det 74fc)/(det A) (for /c = 0,1,..., n). 


C K'nAn 
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It remains to build the polynomial from those coefficients A; = 0,1,..., n. With 
the definition of the vector b we easily obtain 

(First n rows of A) 

0 ... 0 0 ... 0 


(det Ak) k 

Jb • 


Kn 

In the above determinant the entry is in the {k + l)-st column of the matrix. 
So the desired formula 


1 


det A 


^(det Ak)x’^ = 




k=0 


det A 


(First n rows of A) 


1 X x‘^ 


X 


follows from expansion of the matrix by minors along the last row. 


□ 


4. Moment Formulas 

Because our representation for F"(x) effectively depends on the moments we 
now develop a recursive formula for the moments. We do this in an indirect fashion. 
Namely, rather than working with the moments we initially develop a recursion 
formula for the adjusted moments Jlk- We then apply the binomial formula to derive 
the expression for /i^. We choose this indirect approach because it became evident 
during the course of our investigation that the moment calculations for the adjusted 
moments are more concise. (It is possible to do the computations for the exceptional 
moments directly, but such a direct computation turns out to be rather messy.) 

4.1. Adjusted Moments. 

Theorem 4.1. The adjusted moments Jlk = W°'{x) dx 

(a) satisfy the recursion formula 

Jlk+2 = {‘2a + k)Jlk+i + - k)Jlk {keNo), 

(b) and we can start the recursion with 

(4.1) Jlf = e°‘a°‘T{l + a)T{—a, a), and 

= F(a) - 2e"a“F(a + l)F(-a, a), 

where the Gamma function is given by F(x) := t^~^e~^dt and the incomplete 

Gamma function by F(a,a:) := t°‘~^e~^dt for x > 0. 

We simplify notation by writing W for lF“(x). 

Proof. We begin the proof of part (a) by observing two facts. 

First, for functions f,g smooth on [0, oo) the moment functionals satisfy: 

{W',f) = -{W,f) and {gW,f) = {W,fg) 

(here (•, •) denotes the inner product with respect to Lebesgue measure on [0, cxo)). 
Second, we learn from [19] that for a linear operator of the form 

i[y] = a^y” + aiy’ + aoy, 

the related symmetry equation is given by 

022/' + (o's - ai)y = 0, 
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and that it is solved by the weight function (with respect to which the eigen-polynomials 
are orthogonal). 

That is 

a2W' + {a '2 — ai)W = 0 . 

And together with 

(aahh', {x + a)^) = {W, 02 ( 0 ; + af) = - (IT, {a 2 {x + af)') 

= —k (it, a2{x + — (IT, a'2{x -1- a)^) 

we hnd for k eN: 

0 = (oaPT' + (4 - ai)lT, {x + a)^> 

= (aalT', {x + a)’^) + (^IT, {x + a)’^) - (ailT, {x + a)’^) 

(4.2) = -k (IT, 02 ( 2 : + a)^"^) - (ailT, (x + a)^) . 

In our case, with the exceptional Xi-Laguerre expression (1.1) the coefficients are 

2x 


02 = —X and ai = x — a — 1 + 


X + 


or equivalently. 


a 2 = —(x + a) + a and Oi = (x + a) — 2 q! + 1 — 


2a 

X + a 


We substitute these into (4.2), and collect the coefficients of terms of the form 
hfc = (it, (x -f- a)^) for A; = A; — 1, fc, A; -|- 1 to obtain: 

0 = —k (it, (—(x -|- a) + q:)(x + a)^~^) 

'x + a) — 2 q; + 1-^ ) W, (x + a)^ 

X + a J 

= [-ka + 2a]Jlk-i + [k + 2a\ilk - Jlk+i- 
We solve for Jik+i 

Jlk+i = i2a + k- l)pfc q:(2 - k)'flk-i. 

Finally, shifting the index up by one, we see part (a) of the theorem. 

We proceed to prove part (b). We recall the definition of r(x) and r(a, x), and let 

/ OO 

e-^^f-adt, X > 0 

denote the exponential integral function. The two classes of functions are related by 

Ea{x) = x““^r(l — a, x). 

We also have the identity 

(4.3) (a - l ) Ea { x ) = e~^ - xEa - i { x ). 

Our first claim is 


(4.4) 


/•OO g X^l3 

/o (x + a) 


dx = e“ Ei^p{a) r(l + /3), a > 0, fd > —1. 


'OO 
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The claim is established by the following chain of manipulations: 


roo g-x^/3 

/o (x + a) 


dx = J ^\a{t — l))^t ^dt 

pOQ 

= / e“-*(s - afs-^ds 


f*oo poo 


= e 


e ^^{s — aYdtds 


'a J1 
poo poo 


= e^j j e-^'^-^{u/tft-^dudt 

/ oo poo 

e-^^t-^-^dt j e-^u^du. 

As an immediate consequence, we obtain the following expressions for the hrst ad¬ 
justed moment: 

Pi = e“Ei+Q,(a)r(l -h a) 

= r(a)(l -ae“E„(a)) 

= e“a"r(l -|- a)r(—a, a). 

Further, notice that 


(4.5) 




dx= x'^e ^dx = T{a + l). 


□ 


With part (a) for fc = 0, we obtain r(Q; -|- 1) = 2api -|- apo or equivalently 

Po = r(a + l)/a — 2pi = r(a) — 2e"Q;°T(a -|- l)r(— a, a). 

The theorem is proved. 

Equation (4.5) is of interest by itself. 

Remark. R. Milson contributed the evaluation of the moment Jli- 

In the remainder of this subsection, we express the (/c-|-2)-nd adjusted moment pfc +2 
“directly” in two ways. Dehne the matrix 


B„ ■ = 


2a + n a{l — n) 
1 0 


Corollary 4.2. For k >2 


Rk+2 

Rk+1 


r(a + l) mR, 


yn=2 


20-1-1 

1 


Further, we have p 2 = r(Q; -|- 1) and pa = (2a -|- l)r(a -|- 1). 

(The moments po o^nd pi were given explicitly in Theorem 4-1 ■) 
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Remark. The first two moments /xq and Jli occur somewhat disconnected from the other 
moments /x^ for k >2. The forward recursion can be started from k = 2, and that the 
formulas for JIq and Jli contain the incomplete gamma function T. This observation also 
reflects the fact that the exceptional Xi-Laguerre polynomials do not have the degree 
zero polynomial. 


Proof. By the recursion formula in Theorem 4.1 we can set up a discrete dynamical 
system to express 


Pk+2 

= Bu 


Pk+l 




Repeated application yields 


(4.6) 


/^ fc +2 

= 


/^ fc +1 

\t\ ) 

P2 


We recall that p 2 = r(D + 1) by equation (4.5). 
formula in Theorem 4.1 for fc = 1 we have 


Further with the moment recursion 


(4.7) /X 3 — (2q. + l)/i 2 — (2tt + l)r(Q! + 1). 

Substitution into the vector on the right hand side of (4.6) yields the corollary. □ 


We use the standard technique of generating functions (see e.g. [14]) to find an 
explicit expression of the moment. 


Theorem 4.3. The adjusted moments Hk = /o°°(^ + dx are given in hyper¬ 

geometric notation by 

K +2 = (-l)^r(Q; + l)(-a - k)k iFi{-k, -a - k] a) {k E Nq), 
where we used the Pochhammer symbol 


j 1 for n = 0 

\ x{x + 1 ) • ... • (x + n — 1 ) for n > 0. 


(Again, the moments /Xq and /x" cannot be obtained in this fashion, but their values 
are given in Theorem f.l.) 

We check that indeed fTf = r(Q! + l) and JTf = (2a + l)r(Q! + l). The hypergeometric 
lay may consult equation (4.16) below to find an alternative expression for the adjusted 
moments without hypergeometric notation. 


Proof. We begin by re-writing the recurrence relation from Theorem 4.1: 


(4.8) 

to obtain 

/Xfc +2 = (2a fc)/Xfc+i a(l - fc)/Xfc 

{k E No) 

(4.9) 

where 

{k l)uk+i = (2a /c -h l)zxfc - auk-i 

{k E N), 

(4.10) 

fik+2 = kluk {k E No). 
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To see this we first replace := /ifc +2 in (4.8) to see 

Tk = {2a + k)Tk-i + a(l - k)Tk -2 {k = 2, 3,4,...). 

Shifting the index k ^ k + 1 yields 

Tk+i = {2a + k + l)Tk - akTk-i {keN). 

We dehne the seqnence by Tk = k\uk and so 

{k + l)\vk+i = {2a + k + l)k\vk — akluk-i {k E N). 

We divide the equation by k\ to arrive at the desired relation (4.9). 

Our next goal is to write a related hrst order differential equation. In order to achieve 
this, we multiply relation (4.9) by the factor and sum up for fc G N: 

OO OO OO 

(4.11) 'y ^ {k + V)Vk+it^ = ^ ^ (2ci + /c + l)i'kt^ — ^ ^ auk-it^■ 

k=l k=l k=l 

In what follows, t is treated as the independent variable. We dehne the generating 
function 

OO 

(4.12) G{t) :=Y,^kt^- 

/c=0 

Note that the moments can be expressed by 

(4.13) rM = kW = G ^'‘\ 0 ). 

In order to write (4.11) using G(t), we substitute 

OO 

+ l)l^A:+l^^ = G'{t) — Vi 

k=l 

on the left hand side, and 

OO OO 

{2a + 1) ^ ^ kukt’^ = {2a + l)[G{t) — z/q] + tG'{t), as well as 

k=l k=l 

OO OO 

— ^ auk-it’^ = —ctt ^ h’kt'^ = —atG{t). 

k=l k=0 

Apriori, we expect a hrst order inhomogeneous diherential equation. However, when 
we collect the terms without the generating function (i.e. the inhomogeneity) we see 

fi - {2a + l)z/o = P 3 - {2a + l)p 2 = 0 

by Equation (4.7). We obtain the hrst order homogeneous diherential equation 

(4.14) (1 - t)G'{t) + {at -2a- l)G{t) = 0. 

The boundary condition 

(4.15) G'(O) = z^o = h2 = r(a + 1) 
follows from (4.12) for f = 0. 
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It is not hard to see that (near t = 0) the generating function 

G{t) = T{a + 

solves this boundary value problem. 

With the Leibniz rule and elementary simplihcations we derive 


(4.16) «+2 = r(a+l) 


k—m 


m=0 


k 


m 


a^{-a - k+ m)k-m {keNo). 


The theorem now follows from the standard hypergeometric notation by using the 
identity 

1 / I 1 {-a-k)k 

{-a -k + m)k-m = 7 -• 

(We note that the numerator in the latter expression is independent of the summation 
variable m.) □ 

Remark. This method along with a first version of the first order differential equation 
was mentioned to us by M.E.H. Ismail. 

4.2. Exceptional Moments. Starting from the recursion formula in Theorem 4.1 for 
the adjusted moments, we now derive a recursion formula for the exceptional moments 
by using the binomial formula. 

Theorem 4.4. The moments fik = Jq 
(a) satisfy the recursion formula 

k 


Rk+2 — 


m=0 


(2Q! + k) 


k + 1 


m 


— a 


k + 2 


m 


+ {l-k) 


™/im + (1 ~ D)/c/ifc+l 


for k eNq. Specifically, when k = t), we find that p 2 = «(« + l)ho- 
(b) with 

= r(a) - 2e"a“r(a + l)r(-o, a) 

/i“ = —r(Q; + 1) + [2a + l]e“a“r(a + l)r(—a, a). 

Proof. To prove part (b) note that 


p“= / x^W^{x)dx= / {x + afW^{x)dx = fL‘^ 

Jo Jo 

and the first statement follows. To see the second statement 


Jlf= I {x + a)W^{x)dx= I x^W‘^{x)dx + a I x°W^{x)dx = nf + ap.'^. 


So we have 

Substituting (4.1) as well as the formula for /Iq we obtain part (b). 


, .O' 7 TQ ^ . .O 
- /^1 . 
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Part (a) follows from lengthy but elementary computations, which we present for 
the convenience of the reader. The binomial theorem yields 


(4.17) 


= 


1 °° (e 

k 

E 


m=0 

k 


k\ 

I x'^W^(x)dx 

m 


/ifc — ^ 


m=0 


m 


fXrn ^ 


k—ra 


By splitting off the last term see that 


k+2 


/^fc+2 - ^ 


m=0 


k + 2 


m 


fc+i 


lim 


k-\-2—m _ 


E / /c + 2 1 

( _ ) +771 T +k + 2i 


m=0 


m 


that is, 


+k+2 — — ^ ) +7n0t^^‘^ + /ffc+2- 

m=0 


Replacing +k +2 by the recursion formula for the adjusted moments in Theorem 4.1 
and using (4.17) twice, we have 


(4.18) 

(4.19) 


k+1 


+k +2 — ^ _ ) +771 ™ + {2a + k)fik+i + ci(l ~ k)^k 


771=0 

k+1 


m 


I, — — I k + 2 ^ k+2-77l 

+ k + 2 / J \ _ ) +771 Q) 


771=0 


m 


k+1 


+ (2q! + /c) ^ ^ 


771=0 

k 


k + 1 


m 


+771 Cl 


k + 1—771 




m=0 


m 


fJ^Yn ^ 


k—m 


(4.20) 



14 


CONSTANZE LIAW AND JOHN OSBORN 


Finally, we combine the sums. This is done by using 


=a 


k+l 

E 

m=0 
fc+1 

-E 

m=0 


k + 2 


m 


k+l 


flm ^ 


k-\-2—m 


+ 2a^ 


k + 2 


m 


m=0 

fc +1 




fc+1—m 




k + l 


m 


k + l 


+m Cl 


fc+1—m 


=a 


k+l 

E 

m=0 


k + 1\ fk + 2 


m 


m 


+m D 


m=0 
fc+1—m 


m 


+m D 


fc+1—m 


kaj+k+i + D ^ ^ 


m=0 


A; + 1\ fk + 2 


m 


m 


+m D 


fc+1—m 


where the term in square brackets for m = k + 1 evaluated as follows 

'k + 1\ fk + 2^ 


k + l) VA; + 1 


= 2 - (A: + 2) = -A: 


as well as 


fc +1 

^E 

m=0 


k + l 


m 


+m D 


fc+1—m 


k 


m 


m=0 ^ 


-kfik+1 + 


+ «(1 -k)Yl 

m=0 

+ (1-*:)E 

k 


fj^m ^ 


k 


m 


k—m 


k(^^^)+{i-k) 

\ m J \m 

m=0 L \ / \ 

Substitution into (4.18) through (4.20) yields 

fc 


fJjm ^ 


m=0 
fc+1—m 


flm ^ 


fc+1—m 


/^fc +2 — 


— kafik+i + D 


m=0 


/c + 1 \ fk + 2 


m 


m 


fJjm ^ 


fc+1—m 


(4.21) 


+ kfik+1+ 


m=0 


,d+ + (l_;)('= 

m J \m 


flm ^ 


fc+1—m 


and we can easily verify the coefficients in the theorem. 


□ 


5. Alternative Representation of L" with Adjusted Moments 

In Section 4 we saw that the adjusted moments can be computed much more eas¬ 
ily. Here we obtain the more elegant representation for the exceptional Ai-Laguerre 
polynomials. 

We begin by expressing the Ai-Laguerre orthogonal polynomials in terms of powers 
of {x + a). Consider 

n 

Ln{x) = ^ ank{x + aY- 
fc =0 
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As in Section 3, the idea is to express the coefficients ank, A; = 0,1,..., n, via a system 
of n + 1 linear equations Aa = b, where 


a = 

^nO 

e 

and 

b = 

0 


_ ^nn . 




O 


with Kn 7 ^ 0 and where A is given in the next theorem: 


Theorem 5.1. The exceptional Xi-Laguerre polynomials admit the representation 


K{x) = 


det A 


fdet Akj {x + a)^ = 


Kr, 


fc =0 

for n eN, where where 


det A 


(^First n rows of the matrix A 
1 {x + a) (x + ... {x + a) 


A = 


-1 1 0 
Mo F p-i Ml F M2 ■ ■ 
M2 Mo 


0 

Mn F Mn+1 
Mn+2 


/^n+1 


/^2n 


the exceptional moments are given by pk = J^{xFa)’‘W^{x)dx, and where the matrix 
Ak is obtained from A by replacing the {k F l)-st column with the vector b; as is done 
in Cramer’s rule. (In Section 4 we found a recursion formula for the moments fik-) 


Proof. The idea of establishing matrix A is the same way as was for matrix A in 
Theorem 3.1. Since the moments are adjusted to better assimilate Vk{x) = (x + a)^, 
/c > 2, from equations (1.5), the binomial formula is not required and A turns out 
simpler that A. 

The hrst row of A also simplihes: As before, we use the exceptional condition (2.1). 
We obtain 

n 

(L“)' (x) = ^ kank{x F a)^~^ 

k=l 


and so 


-^n(-«) = Ono as well as (L“)' (-a) = a^i. 


We have 

0 = (L")' (-a) - T"(-a) = Oni - a„o. 

The entries in the first row of A follow. 

The other n conditions are, again, obtained via orthogonality. So we have 


('^n ; ^ndnk {k 1, . . . , 77.). 
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For A; = 1, we observe as before 


, k=0 


ank{x + a)^ , X + a + 1 

I {x + a + l)W°‘{x) dx 


Kk=0 


^ ^ ^nki^f^k /^/c+l) • 


k=0 


So the second row of the matrix A is the first vector in the dot product 

[ /io + /^1 /^1 + /^2 • • • /^n + f^n+1 ] ■ a = Kndnl- 

When n > 2 we consider 2 < s < n. Again we compnte 

K-ndns — i Vs)^ra 

ank{x + a)^ , (x + ay 

I k =0 / lyo 

n 

^ ^ (^nk f^k+s- 
k =0 

The corresponding coefficients fill the (s + l)-st row of the matrix A. 
Again, the matrix A is invertible and we apply Cramer’s rnle 

(for fc = 0,1,..., n). 

det Ak 


ttnk = [det Akj / [det A 
With the definition of the vector b: 


(First n rows of A) 

0 ... 0 (x + a)^ 0 ... 0 




(x + a) , 


and the desired formnla follows from expansion of the matrix by minors along the last 
row. □ 

Remark 5.2. In order to embed this representation into the literatnre, we relate to the 
normalization used in [6] and [17]: The choice 

Kn= (-!)"(«+ n)F(a + n-l) 

yields the same normalization as in [6] and [17]. 

Indeed, there the leading coefficient of T"(x) is given by (—1)"'/(?7, — 1)! and with 
this normalization they obtained 

2 F(q; + n - l)(a + n) 


|L“IF = 


(n — 1)! 
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And so by the Gram-Schmidt orthogonalization we have 

= (L“(x), {x + a)”) = (n- l)!/(-l)”||L“||2 = (-l)’"(a + n)r(a + n - 1). 

This relation was pointed out to us by R. Milson. 

Finally, we verify for n = 1 and n = 2 that the polynomials in Theorem 5.1 indeed 
agree with (1.2) and (1.3), respectively, as well as the normalization claimed in the 
latter remark. 

Example. For n = 1 compute 


-1 1 
1 X + a 


= —X — a — 1, 


which is a scalar multiple of (1.2). And, moreover, normalizing the formula in Theorem 
5.1 according to Remark 5.2 we obtain 

— {a + l)F(a) 


det A 


-1 1 
1 X + a 


-{—X — a — 1) = —X — a — 1, 


{—a — 1)F(«) 

the leading term of which is in agreement with the leading term in Remark 5.2, since 
for n = 1 we have (—l)"'/(n — 1)! = —1. 

Example. Take n = 2. By co-factor expansion along the last row we evaluate 

-1 1 0 
— ho + hi hi + h2 h2 + h3 

1 X + a (x -h a)^ 

= h2 + hs T (^ T ci)(h2 + hs) + (31 -|- cy)'^[—jj.Q — 2/ii — /i2]. 

For A; = 0 in part (a) of Theorem 4.1 the recursion relation reduces to 

P 2 — 2,a.p.i -|- ci/Uo, 

so that 

—ho ~ 2/ii = —h 2 /<T- 

And recall that the recursion relation for k = 1 reduces to 

yhs = (2a -h l)h 2 . 

With this we have up to normalization 

^2 (x) = (2a -|- 2 )yU 2 -|- (x -1- a)(2a -|- 2 )yU 2 — (x -|- a)^(l/a -|- l)h 2 - 
Factoring out 


1 a-|-l.^ 

- + 1] 1^2 = -h2 

a a 


we obtain 


= 


H“ 1 
a 

“f" 1 
a 

Q H“ 1 
a 


Ji 2 [—2a — 2a(x -|- a) -|- (x -|- a)^] 
yU 2 [—2a — 2ax — 2a^ -1- x^ -|- 2ax -|- a^] 
Ji 2 [x^ — a^ — 2a] . 
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This is in agreement with the degree polynomial given by (1-3), since they both span 
the same eigenspace. 

Again, we compare this to the Remark 5.2 about normalization. It is tedious but 
elementary to compute 

K 2 _ —{a + 2) 

detA + 3 + 2 / a) 

for K 2 as in Remark 5.2. Now, the leading term of (with the normalization 

described in Theorem 5.1 and with K 2 from Remark 5.2) is given by 

/ a + \ / Kn ^ + 1)(® + 2) _ 1 

V ^ ^ \ det AJ + 3a + 2 

And, again, this is in agreement with the predicted leading term coefficient in Remark 
5.2, since for n = 2 we have (—— 1)! = 1. 

Acknowledgement. We would like to thank K. Busse and L.L. Littlejohn for useful 
discussions. 
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